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Abstract 

We analyze the phenomenon of system size stochastic resonance in a simple spatially extended 
system by exploiting the knowledge of the nonequilibrium potential. We show that through the 
analysis of that potential, and particularly its "symmetry", we can obtain a clear physical in- 
terpretation of this phenomenon in a wide class of extended systems, and also analyze, for the 
same simple model, the effect of a general class of boundary conditions (albedo) on this kind of 
phenomena. 
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Among the many noise induced phenomena extensively studied during the last few 
decades like: stochastic resonance jl|, noise- induced transitions Q|, noise- induced phase 
transitions , noise-induced transport jj, [J , noise- induced limit cycles ; stochastic res- 
onance (SR) detaches as one that attracted -and still attracts- the attention of a large 
number of researchers, probably due to its interest from a technological as well as a biologi- 
cal points of view. There is a broad range of phenomena for which this mechanism can offer 
an explanation as has been put in evidence by many reviews and conference proceedings [lj]. 

The fingerprint of the SR phenomenon is the enhancement of the output signal-to-noise 
ratio (SNR) caused by the injection of an optimal amount of noise into a periodically -or 
even aperiodically- driven nonlinear system. Such enhancement is the result of a cooperative 
effect arising from the interplay between deterministic and random dynamics in a nonlinear 
system. In almost all the studies of SR, the relevant control variable of the phenomenon 
was the noise intensity, while the system's size didn't play any relevant role. However, some 
recent studies on biological models of the Ho dgkin- Huxley type 0, 0] have shown that ion 
concentrations along biological cell membranes presents intrinsic SR-like phenomena as the 
number of ion channels is varied. A related result |9] shows that even in the absence of 
external forcing, the regularity of the collective firing of a set of coupled excitable FitzHugh- 
Nagumo units results optimal for a given value of the number of elements. From a physical 
system point of view, the same phenomenon -that has been called system size stochastic 
resonance (SSSR)- has also been found in an Ising model as well as in a set of globally 
coupled units described by a 4 theory . It was even shown to arise in opinion formation 



models 



0. 



In a recent series of papers |12J, |13|, |14|, [la, |16( the phenomenon of SR in extended systems 



was studied exploiting the concept of nonequilibrium potential (NEP). This potential is a 
special Lyapunov functional of the associated deterministic system which, for nonequilibrium 
systems, plays a role similar to that played by a thermodynamic potential in equilibrium 
thermodynamics Q|. Such a NEP, closely related to the solution of the time independent 
Fokker-Planck equation of the system, characterizes the global properties of the dynamics, 
that is: attractors, relative (or nonlinear) stability of these attractors, height of the barriers 
separating attraction basins, allowing to evaluate the transition rates among the different 



attractors. In 



I'H llflj - this approach was applied to the global stability analysis of some 



react ion- diffusion systems. A kind of "mini-review" on the application of this approach to 
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he study of SR in one, two, and three species reaction-diffusion systems could be found in 
16^ . However, in spite of its potentiality, this kind of approach has not been exploited for 
the study of SSSR. 

In this paper we present an analysis of the SSSR phenomenon in a simple spatially 
extended system, exploiting previous results obtained using the notion of NEP within the 
context of a simple reaction-diffusion model. The specific model we shall focus on, with a 
known form of the Lyapunov function, corresponds to a one-dimensional, one-component 



model 
models 



20L2 



2jJ that, with a piecewise linear form, mimics general bistable reaction-diffusion 
In particular we will exploit some ^of the results for the influence of general 

as well as previous studies of the NEP 



boundary conditions (called "albedo") found in 
3,Q and of SR [l2, Q. 

The particular non-dimensional form of the model that we work with is |l2 



13. 
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dtc 



(1) 



We consider here a class of stationary structures (p(y) in the bounded domain y G [— Vl^Vl] 



with albedo boundary conditions at both ends, ^ 



v=±vl 



0(±?Al), where k > is the 



albedo parameter. These are the spatially symmetric solutions to Eq.(^) already studied in 



22 . The explicit forms of these stationary patterns -not shown here- are given by Eq. (9) 
of 3 (see also Fig. 4 in 22 1). 



The double- valued coordinate y c , at which 



is given by 



He 
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z*y(k, y L ) ± a/ z 2 ^(k, y L ) 2 + 1 - ¥ 
1 + k 



(2) 



with 7(/c,y) = sinh(y) + k cosh(y), and z — 1 — 2(f> c /(ph (—1 < z < 1). 

When exists and yf < yL, this pair of solutions represents a structure with a central 
"excited" zone {(f) > <f) c ) and two lateral "resting" regions (<f) < <p c ). For each parameter set, 
there are two stationary solutions, given by the two values of y c . Figure 5 in 22\ depicts the 
curves corresponding to the relation y c /yL vs. k, for various values of 4> c /(f)h- 

It has been shown (0 that the structure with the smallest "excited" region (with y c = y+, 
denoted by (fi u (y)) is unstable, whereas the other one (with y c = y~ ', denoted by 4>i(y)) is 
linearly stable. The trivial homogeneous solution <fio(y) = (denoted by 0o) exists for 
any parameter set and is always linearly stable. These two linearly stable solutions are 
the only stable stationary structures under the given albedo boundary conditions. We will 



concentrate on the region of values of z, jjl and k, where both nonhomogeneous structures 
exist. 

For the system with the albedo b.c. that we are considering here, the NEP reads 



r[</>,k, yi ]= I I - I [-</>' + <t> h e{cj>' -<j> c )\ d<P' + -(d y <P(y,t)) 2 \dy+ -<P(y,tf 
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(3) 



Replacing the explicit forms of the stationary nonhomogeneous solutions (see for instance 
Eq.(9) in we obtain the explicit expression [12I llflj] 

J* = T[K,x, k, y L ) = -eg y±z + eg sinh^) ~f\ (4) 

while for the homogeneous trivial solution <p Q = 0, we have instead T[4>q, k,yi] = J 70 = 0. 

Figure 1, part (a) depicts k, y{\ as a function of the system size yL, for a fixed 
albedo parameter k, and a fixed value of the ratio (pc/^h (i-e. fixed value of z). The curves 
correspond to the nonhomogeneous structures, JF ± , whereas the horizontal line stands for 
JF°, the NEP of the trivial solution. We have focused on the bistable zone, the upper branch 
being the NEP of the unstable structure, where T attains a maximum, while in the lower 
branch (for = O or — 0i)> the NEP has local minima. We see that when yi becomes 
small, the difference between the NEP for the states (f> u {y) and <f>i(y) reduces until, for 
yL ~ 0.72, they coalesce and, for even lower values of yL, disappear. For completeness and 
latter use, in part (b) of Fig. 1 we show J-[(f>, k,yi\ but now as a function of k, for a fixed 
value of yL and the same value of z. Here we see that the initial large difference between 
the NEP for the states </> u (y) and 4>i(y) reduces for increasing k until, for k — > 00, the values 
for Dirichlet b.c. are asymptotically reached. 

It is important to note that, since the NEP for the unstable solution <fi u is always positive 
and, for the stable nonhomogeneous structure <f>i, T < for yL large enough, and T > for 
small values of yL, the NEP for this structure vanishes for an intermediate value yL = y\ °f 
the system size. At that point, the stable nonhomogeneous structure <f>x(y) and the trivial 
solution 4>o(y) exchange their relative stability. 

In order to account for the effect of fluctuations, we include in the time-evolution equation 
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of our model (Eq.((TJ)) a fluctuation term, that we model as an additive noise source 
yielding a stochastic partial differential equation for the random field <j)(y, t) 

d t cj>{y, t) = d 2 yy <p - <p + 0^(0 - C ) + ^(y, *)■ (5) 
4 



FIG. 1: NEP evaluated at the stationary solutions <j)o(y), </>i(y) and 4> u {y)- Part (a) T\<\>, k,y£\ vs. 
Ul, with k = 3. Part (b) .F[</>, k,yi\ vs. k, with yi = 1.2. In both cases ^c/^h = 0.193. 



We make the simplest assumptions about the fluctuation term £(y,t), i.e. that it is a 
Gaussian white noise with zero mean and a correlation function given by: (£,{y, t) £,(y', t')) = 
2 7 S(t — t') 8(y — y'), where 7 denotes the noise strength. 



As was discussed in 
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13, 



14. 
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161 ] . using known results for activation processes in 



multidimensional systems |2J] , we can estimate the activation rate according to the following 
Kramers' like result for (r), the first-passage-time for the transitions between attractors, 

~AJ*[<j>, k] 



(Tj) = r exp 
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(6) 



where A.P[0, k, y{\ = F[4> u (y), k, yi] - F[(/>i(y), k, yi] (i = 0,1). The pre-factor r is usually 
determined by the curvature of ^[(f), k, y{\ at its extreme (minima) and typically is, in 
one hand, several orders of magnitude smaller than the average time (r), while on the 
other does not change significatively when changing the system's parameters. Hence, in 
order to simplify the analysis, we assume here that r is constant and scale it out of our 



results. The behavior of (r) as a function of the different parameters (k, <p c ) was shown in 

We assume now that the system is subject (adiabaticallvV to an external harmonic vari- 
ation of the parameter < f> r : <t> r (t) = <j> r + 8(j) c cos{uit) [3, 13, and exploit the "two-state 



approximation" lj as in 
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16| . For all details on the general two-state approxima- 
tion we refer to 

Up to first-order in the amplitude 5(p c (assumed to be small in order to have a sub 
threshold periodic input) the transition rates Wj take the form 

A^[(f>,k,y L y 



Wi 



r exp 
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(7) 



where 

Ar[4>,k,y L ] = AJ*[& k, y L ] + S^J 9 ^^ ^ Vl] ) cos(arf). (8) 

V 0<p c / <t>c=<j)* c 

This yields for the transition probabilities 

1 / 5(p c \ 

Wi^-\jii^oti—coa{ut)J, (9) 

where ^ « exp(-A.P[0, fc, y L ]) and a { « ±/i<^|f U„ (i = 1,2). Using Eq. ®, |# can 
be obtained analytically. 

These results allows us to calculate the autocorrelation function, the power spectrum and 
finally the SNR, that we indicate by R. The details of the calculation were shown in ^| 
and will not be repeated here. For R, and up to the relevant (second) order in the signal 
amplitude b~<fi C) we obtain 1^ 

R= ?r (a 2 ^ + ai fi 2 ) 2 

We have now all the elements required to analyze the problem of SSSR. 

Figure 2 shows the typical behavior of SR, but now -in the horizontal axis- the noise 
intensity is replaced by the the system length y^, for fixed values of k, 7 (the noise intensity) 
and c /0h- Such a response is the expected one for a system exhibiting SSSR. Within the 
context of NEP, it results clear that, in this kind of systems, the phenomenon arises due to 
the breaking of the NEP's potential symmetry. That is: both attractors change their relative 
stability due to the variation of yi as shown in Fig. l.a. Hence, within this framework, SSSR 
arises as a particular case of the general discussion in jis| . 

Let us now change the point of view. In Fig. 3 we show the curves of the SNR as a 
function of k, while keeping fixed values of yi, and z. When k is not too large, indicating a 
high degree of reflectiveness at the boundary (or a reduced exchange with the environment), 
we see that the SNR changes for k varying from low to larger values. Remember that a large 
value of k indicates that the system boundaries become absorbent. Also, the robustness of 
the systems' response when changing k, a parameter that somehow indicates the degree of 
coupling with the environment, is apparent. According to the previous argument -about 
the breaking of NEP's symmetry- from Fig. l.b this is the expected result. 

As a final remark, let us consider one of the models discussed in from the point of 
view of the above indicated approach. The model we refer to is described by a set of coupled 
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FIG. 2: SNR, vs. y L , for k = 3., 7 = 0.1 and </> c/ % = 0.193. 
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FIG. 3: SNR, vs. k for y L = 1.2, 7 = 0.1 and <f> c /<f> h = 0.193. 
nonlinear bistable oscillators 

N 

Xj = Xj -tf + jj J2( Xk ~ x ^> + V^7^j(*) + /i(*)> 
fc=i 

O Xj 

with = Acos(cjt), {x} = (x 1 , x 2 , .., Xjv), and 

AT AT 

C/({x},t) = ~ Acos(o;t) 



Xj 



Xj 



N 

= Uq({x}) - Acos(ut) 

3=1 

N , 4 2\ iV JV N 

3=1 ^ ' 3=1 k=l 3=1 

Due to the structure of Eq. it is clear that the potential function in Eq. (|T2j) is the NEP 
of this problem and -for A = 0- that the stationary distribution of the multidimensional 
Fokker-Planck equation associated to Eq. (fTTjl is 

PU{x}) « exp V (is) 
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Clearly, this potential has two attractors corresponding to x\ = x^ = ■■■ = xn = ±1, and a 
barrier separating them at X\ = x 2 = ■■■ = x^ = 0. 

Now, exploiting the same scheme as before, but reduced to the symmetric case (as both 
attractors have the same "energy" ) , we get 

SNR w exp ( " JJ J w — exp ( v 1 \ , (14) 



where X is a "collective coordinate" (that can be approximately interpreted as X ~ 
jjYl!j=i x j)-> an d A-uo(^) = uq{x = ±1) — Un ix = 0). This SNR clearly shows similar 
SSSR characteristics as those described in . However, as in this situation the NEP's 
symmetry is retained when varying N, and we could speak of a genuine SSSR. 

The above results clearly show that the "nonequilibrium potential" , even if not known 
in detail offers a very adequate framework to analyze a wide spectrum of noise in- 
duced phenomena in spatially extended or coupled systems. Within this framework the 
phenomenon of SSSR looks, as in other aspects of SR in extended systems , clS db nat- 
ural consequence of the breaking of the symmetry of the NEP. In addition, we have seen 
that through the variation of its coupling with the surroundings, a system can increase the 
robustness of its response to an external signal, opening new possibilities of analyzing and 
interpreting the behavior of biological systems. 

An important conclusion to be drawn from the identification of the Lyapunov functional 
with a "thermodynamical-like potential", is that for a wide range of parameters where 
the system is essentially bistable (with both attractors not necessarily having the same 
"energy"), the problem admits a one-dimensional analogue 



HQ 



16. 18. 191. This feature is in 



contrast with the infinite dimensional character of the whole function space, and has been 



exploited to strongly simplify the analysis. 
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